(N 

o 

(N 



Employing pre-stress to generate finite cloaks for antiplane elastic waves 
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It is shown that nonhncar clastic pre-stress of nco-Hookcan hypcrclastic materials can be used as a mechanism 
to generate finite cloaks and thus render objects near-invisible to incoming antiplane elastic waves. This 
approach appears to negate the requirement for special cloaking metamaterials with inhomogeneous and 
anisotropic material properties in this case. These properties are induced naturally by virtue of the pre-stress. 
This appears to provide a mechanism for broadband cloaking since dispersive effects due to metamaterial 
microstructure will not arise. 

PACS numbers: 46.40.Cd, 62.20.F-, 62.30. -t-d 
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Interest in cloaking theory and its practical realiza- 
tion has grown significantly since the early theoretical 
work of Leonhardt^ and Pendry et al.^ in optics and elec- 
tromagnetism respectively. Methods have been largely 
based on the idea of coordinate transformations^ which 
motivate the design of cloaking metamaterials that are 
able to guide waves around a specific region of space. 
Since the early work, research has focused on the possi- 
bility of cloaking in the contexts of acoustics^"^ , surface 
waves in fiuids^, heat transfer^, fluid flow^ and linear 
elastodynamics^°"^'^ and it is the latter application which 
is the concern of the present letter. In particular it was 
shown in^" that elastodynamic cloaking is made difficult 
due to the lack of invariance of Navier's equations un- 
der general coordinate transformations which retain the 
symmetries of the elastic modulus tensor. A special case 
is that of flexural waves in thin plates^*^. Invariance of 
the governing equations can be achieved for a more spe- 
cific class of transformations if assumptions are relaxed 
on the minor symmetries of the elastic modulus tensor 
as was described for the in-plane problem in^^. Cosserat 
materials were exploited in^'^. 

As noted in^^, another special case for elastodynam- 
ics is the antiplane elastic wave problem, where cloaking 
can readily be achieved from a cylindrical region (using 
a cylindrical cloak) in two dimensions by virtue of the 
duality between antiplane waves and acoustics in this 
dimension. Consider an unbounded homogeneous elas- 
tic material with shear modulus fiQ and density po and 
introduce a Cartesian coordinate system {X, Y, Z) and 
cylindrical polar coordinate system {R, Q, Z) with some 
common origin O. Planar variables are related in the 
usual manner, X = R cos Q,Y = RsinQ. Suppose that 
there is a time-harmonic line source, polarized in the Z 
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direction and located at (i?o, ©o); with circular frequency 
a; and amplitude C (which is a force per unit length in 
the Z direction). This generates antiplane elastic waves 
with the only non-zero displacement component in the Z 
direction of the form U = ^W{X,Y)cxp{~iujt)]. The 
displacement W is governed by 

Vx-(moVxVF)+pow2^ = C(5(X-Xo), (1) 

where Vx is the gradient operation in the "untrans- 
formcd" frame, X = {X,Y) and Xq = {Xo,Yo). 

The assumed mapping for a cloak for antiplane waves 
(cf. acoustics) expressed in plane cylindrical polar coor- 
dinates, takes the form 



5(i?), 



e. 



Z, for < i? < i?2 



(2) 



and the identity mapping for all R > R2 for some chosen 
monotonically increasing function g{R) with g(0) = ri G 
[0,i?2], g{R2) = i?2 e K such that i?2 < ^0, i-e. the line 
source remains outside the cloaking region. The cloaking 
region is thus defined by r S [?'i,?'2] where r2 = i?2- We 
use upper and lower case variables for the untransformcd 
and transformed problems respectively. Under this map- 
ping the form of the governing equation (1) remains un- 
changed for i? = r > i?2, whereas for < i? < i?2, 
corresponding to the transformed domain ri < r < R2, 
the transformed equation takes the form (in transformed 
cylindrical polar coordinates r,9 — Q) 
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+ d{r)uj^w = (3) 



where (see eqs. (26), (27) in^) 
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Hence, both the shear modulus and density must be inho- 
mogeneous and the shear modulus must be anisotropic. 



Material properties of this form cannot be eonstructed 
exactly since the shear modulus fj,g becomes unbounded 
as r — ^ ri (the inner boundary of the cloak). In this limit 
the density behaves as d = {pcri)~^ poR^~P + . . . where 
p, c > define the mapping in the vicinity of the inner 
boundary according to r = ri + cRf + . . . as i? — s- 0. 
In practice of course approximations are required as de- 
scribed in c.gJ'^^'^^. Note that, as expected^'^, the total 
mass is conserved since, regardless of the mapping, the 
integral of the density d{r) over r £ [ri,r2] is irR'^po- 

In^^ a new method to generate elastic cloaks was pro- 
posed which used the notion of nonlinear pre-stress. That 
this was possible was due to the fact that the antiplane 
wave field scattered from a cylindrical cavity is invari- 
ant under pre-stress for an incompressible neo-Hookean 
material. Scattering coefficients in the deformed config- 
uration depend only on the initial cavity radius i?i and 
therefore provided that this is small compared with the 
incident wavelength, scattering from the infiated cavity 
of radius ri will be negligible regardless of the relative 
size of ri and the incident wavelength. Therefore we 
can conclude that an object placed inside the infiated 
cavity region would be near-invisible (i.e. cloaked) upon 
choosing Ri appropriately. In^^ the pre-stress affected 
the entire elastic domain however and therefore its influ- 
ence was felt by both the source and receiver. In this 
letter we show how this theory may be adapted in order 
to create a finite cloak by means of an axial stretch. 

With reference to Fig. 1, let us consider an elastic ma- 
terial within which is located a cylindrical cavity of radius 
i?2- Let us assume that the density of this medium is po 
and its axial shear modulus (corresponding to shearing on 
planes parallel to the axis of the cylindrical cavity) is /io- 
Additionally we take a cylindrical annulus of isotropic in- 
compressible neo-Hookean material with associated shear 
modulus /x and density p and with inner and outer radii 
Ri and i?2 respectively with Ri <€. R2- The exact nature 
of this latter relationship will be described shortly. We 
shall consider deformations of the cylindrical annulus in 
order that it can act as an elastodynamic cloak to incom- 
ing antiplane elastic waves. We deform the material so 
that its inner radius is significantly increased (to ri) but 
its outer radius R2 remains unchanged. The deformed 
cylindrical annulus can then slot into the existing cylin- 
drical cavity region within the unbounded (unstressed) 
domain. We choose p and p so that subsequent waves 
satisfy the necessary continuity conditions on ?' = i?2. 

The constitutive behaviour of an incompressible neo- 
Hookean material is described by the strain energy 
function^® 



W^^{\l + \l + Xl-'i) 



(5) 



where A,, j 



J, J — r,9, z are the radial, azimuthal and ax- 
ial principal stretches of the large deformation to ensue. 
We consider the initial deformation of the cylindrical an- 
nulus domain as depicted in Fig. 1. Since the material 



FIG. 1. The incompressible neo-Hookean cylindrical annulus 
is pre-stressed as depicted on the right. This annulus then 
creates a cloak when slotted into a cylindrical cavity in an 
unbounded elastic medium, as illustrated on the left. 



is incompressible, the deformation is induced either by 
applying a uniform axial stretch L or a, radial pressure 
difference Po — Pi where Po and pi denote the pressures 
applied to the outer and inner face of the cylindrical an- 
nulus respectively. The ensuing deformation is described 
via the relations 



R^R{r), 



e 



Z = z/L, 



(6) 



where (R,Q,Z) and (r,6,z) are cylindrical polar coor- 
dinates in the undeformcd and deformed configurations. 
Note the convention introduced in (6), i.e. that upper 
case variables correspond to the undcforined configura- 
tion whilst lower case corresponds to the deformed con- 
figuration. This is analogous to the notation used for 
untransformcd and transformed configurations in (2). 
The principal stretches for this deformation are 

dr 1 r 

^'^ = dE = i?v)' ^'^W)' ^'^^' ^^^ 

For an incoinpressible material XrXgXz = 1, implying 



R{r) = ^L{r^ + M), 



(8) 



where M = R2{L~^ — 1) is a constant determined by 
imposing that the outer wall of the cylindrical annulus 
remains fixed, i.e. i?(i?2) = R2- The deformation (8) 
is easily inverted to obtain r{R). Given incomprcssibil- 
ity and the fixed outer wall of the annulus, in order to 
induce this dcforination we may either (i) prescribe the 
axial stretch L which then determines the deformed in- 
ner radius ri and the radial pressure difference required 
to maintain the deformation or (ii) prescribe the radial 
pressure difference which then determines the deformed 
inner radius ri and the axial stretch L. 

We shall discuss the radial pressure difference shortly 
but either way we can obtain L and thus feed this into (8) . 
Imposing the requirement that R{ri) = Ri and using the 
form of M gives rise to the useful relation 



Rl-rl 



(9) 



The Cauchy stress for an incompressible material is 
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(10) 



where W is the neo-Hookean strain energy function in- 
troduced in (5), F is the deformation gradient, I is the 
identity tensor and Q is the scalar Lagrange multiplier 
associated with the incompressibility constraint. 

Only diagonal components of the Cauchy stress are 
non-zero, being given by (no sum on the indices) 



the deformation (6) and (8). We can obtain the gov- 
erning equation using the theory of small-on-large^®. It 
was shown in^^ that the wave in this region satisfies (3) 
but now where iJ-rir) and ^J.e{r) are defined in (12) and 
d{r) = p, and note that we have made the necessary cor- 
rections in order to include the axial stretch L which was 
not considered in^^. Note in particular that the density 
is homogeneous inside the cloak region. 

Let us introduce the identity mapping for r > R2 and 

R^ = L{r^ + AI), e = e, for n < r < R2 (15) 



T,. 



ioT j = r,9,z, where 
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(12) 



The second and third of the static equations of equi- 
librium div T = (where div signifies the divergence 
operator in the deformed configuration) merely yield 
Q = Q{f)- The remaining equation 



dT„ 
dr 



-{Trr — Tee) — 0, 



(13) 



can be integrated using (11)-(12) to obtain Q{r). Writing 



T . 



I r=i?,2 



Poi ^•rr\^_^^ 



L{Pi -Po) _ 1 
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-pi we find 



Rl 



= :r^ 1 - ^ + log -^ 
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(14) 



Given L and thus ri via (9), this equation prescribes the 
required pressure difference. 

Now assume that the cylindrical annulus has been 
pre-stressed in an appropriate manner and slotted into 
the unbounded elastic material with perfect bonding at 
r = i?2. We consider wave propagation in this medium 
given a time-harmonic antiplane line source located at 
(i?Oj0o) with Rq > i?2. In r > i?2 the antiplane 
wave with corresponding displacement which we shall de- 
note by w{r,9), is again governed by (1). In the region 
Ti < r < i?2, the wave satisfies a different equation since 
this annulus region has been pre-stressed according to 

I 



which corresponds to the actual physical deformation 
(8). Finally define W{R,Q) = w{r{R),e{&)). It is then 
straightforward to show that the equation governing wave 
propagation in the entire domain i? > i?i is (1), provided 
that we choose fi = Ljjlq and p = LpQ. These relations 
ensure that the wavenumbers in the exterior and cloak 
regions are the same and they also maintain continuity 
of traction on R = R.2. Furthermore since (15) corre- 
sponds to the actual deformation, the inner radius ri 
maps back to Ri . Therefore with the appropriate choice 
of cloak material properties, the scattering problem in 
the undeformed and deformed configurations are equiv- 
alent. We can therefore solve the equation in the un- 
deformed configuration and then map back to the de- 
formed configuration to find the physical solution. De- 
composing the solution into incident and scattered parts 
W = Wi +Ws, "^e have Wi = ■^llo{KS) where we 
have defined the wavenumber K via K^ = pi^bj^ j p^ and 
S = ^{X-X^f + iy-Y^f. Here H„ = H^,!) is the 
Hankel function of the first kind of order n. The scat- 
tered field is written in the forni"'^^ 



Ws{R) 






-i)"a„H„(iri?)e 



m(e-eo) 



(16) 



Satisfaction of the traction free boundary condition on 
R = R\ gives a„ . We want the wave field with respect to 
the deformed configuration, so we map back in order to 
find w = Wi + Ws- The incident wave is most conveniently 
determined by using Graf's addition theorem in order to 
distinguish between the regions r < Rq and r > Rq, as 
was described in^^. The incident and scattered fields are 
then, respectively. 
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(18) 
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FIG. 2. Cloaking of antiplane shear waves. Line source is 
located at Kr — KRo — Svr, Oo = 0, shown as a white 
circle. Upper left: A region of (nondimensionalized) ra- 
dius Kri = 2n is cloaked using a classic linear elastic cloak 
g{R) = r-i + R{^^^) in 2tt < Kr < An. Upper right: 
Scattering from a cavity of radius KRi — 27r/20 in an 
unstressed medium. Lower left: A "pre-stress" cloak in 
27r < Kr < An generated from an annulus with initial inner 
radius KRi — 27r/20. Lower right: Scattering from a cavity 
with radius KRi = 27r in an unstressed medium. Scattering 
and the shadow region presence in the latter is significant, 
as compared with that for an equivalent sized cavity for the 
"pre-stress" cloak. 



The key to cloaking is to ensure that the scattered field 
is small compared with the incident field, i.e. a„ <^ 1. 
Note from (18) that a„ are solely dependent on the initial 
annulus inner radius Ri (and source distance Ro) but are 
independent of the deformed inner radius ri . Therefore 
we must choose i?i such that KRi <Si 1 which will ensure 
negligible scattering. We illustrate with some examples 
in Fig. 2, showing that the "pre-stress" cloak appears to 
work well. 

In conclusion, we have shown how a finite cloak for 
antiplane elastic waves can be generated by employing 
nonlinear pre-stress of an incompressible neo-Hookean 
hypcrelastic material. The performance of the cloak is 
limited only by the size of the initial radius of the cylin- 
drical cavity inside the annulus region. The anisotropic, 
inhomogeneous material moduli in the cloaking region, 
defined by (12), are induced naturally by the pre-stress 
and therefore exotic metamaterials are not required. Dis- 
persive effects, which naturally arise in metamaterials 
due to their inherent inhomogeneity at some length scale, 
will not be present in the pre-stress context and we also 



note that the density of the cloak is homogeneous. In 
order to achieve the required pre-stress, a radial pressure 
difference is required across the cylindrical annulus. It 
would be inconvenient to prescribe Po on the outer face. 
However, since we only need a pressure difference we can 
prescribe pi with po = 0, ensuring the prescribed defor- 
mation and eliminating this difficulty. The incompress- 
ible nco-Hookcan model is an approximation to reality, 
holding in general for rubber-like materials and moder- 
ate deformations. If the material is not neo-Hookean, 
invariance of the scattering coefficients is not guaran- 
teed in general and therefore similar exact results will 
not hold. However it would be of interest to ascertain 
whether scattering from inflated cavities in other hyper- 
elastic pre-stressed media is still significantly reduced as 
compared with an equivalent sized cavity in an unstressed 
medium. 

In closing we remark that one of the fundamental ad- 
vantages of the pre-stress approach is that pre-stress gen- 
crates equations with incremental moduli (analogies of 
the elastic moduli) which do not possess the minor sym- 
metries. Therefore this approach can be used for cloaking 
in the more general elastodynamic setting where classical 
linear elastic materials cannot be used^^. 

ACKNOWLEDGMENTS 

The work of ANN was supported by NSF and ONR. 
TS is grateful to EPSRC for funding his PhD research 
studentship. 

lU. Lconhardt, Science 312(5781), 1777-1780 (2006). 

2 J.B. Pendry, D. Schurig & D.R. Smith, Science 312(5781), 1780- 

1782 (2006). 
^A. Greenloaf, M. Lassas & G. Uhlmann, Math. Res. Lett. 10, 

685-693 (2003). 
"S.A. Cummer & D. Schurig, New J. Phys. 9 45, 8 pp (2007). 
^H. Chen & C.T. Chan Appl. Phys. Lett. 91, 183518 (2007). 
•^A.N. Norris, Proc. Roy. Soc. A 464, 2411-2434 (2008). 
^M. Farhat, S. Enoch, S. Guenneau, & A.B. Movchan, Phys. Rev. 

Lett. 101, 134501 (2008). 
«J.Y. Li, Y. Gao & J.P. Huang, J. Appl. Phys. 108, 074504 (2010). 
Sy.A. Urzhumov & D.R. Smith, Phys. Rev. Lett. 107, 074501 

(2011). 
i"G.W. Milton, M. Brianc & J.R. Willis, New J. Phys. 8, 248-267 

(2006). 
i^M. Brun, S. Guenneau & A.B. Movchan, Appl. Phys. Lett. 94, 

061903 (2009). 
^^A.V. Amirkhizi, A. Tehranian & S. Ncmat-Nasser, Wave Motion 

47, 519-536 (2010). 
13 A.N. Norris & A.L. Shuvalov, Wave Motion 48, 525-538 (2011). 
^''M. Farhat, S. Guenneau, S. Enoch & A.B. Movchan, Phys. Rev. 

B 79, 033102 (2009). 
l^D. Schurig, J.J. Mock, B.J. Justice, S.A. Cummer, J.B. Pendry, 

A.F. Starr & D.R. Smith, Science 314(5801), 977-980 (2006). 
if'S. Zhang, C. Xia, & N. Fang, Phys. Rev. Lett. 106, 024301 

(2011). 
I'^W.J. Parnell, Proc. Roy. Soc. A 468, 563-580 (2012). 
i*R.W. Ogden, Nonlinear elasticity. Dover, New York (1997). 
i^A.N. Norris & W.J. Parnell, Submitted to Proc. Roy. Soc. A. 



